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ABSTRACT 

We report the detection and measurement of the absolute brightness and spatial fluctuations of 
the cosmic infrared background (CIB) with the AKARI satellite. We have carried out observations 
at 65, 90, 140 and 160 /mi as a cosmological survey in AKARI Deep Field South (ADF-S), which 
is one of the lowest cirrus regions with contiguous area on the sky. After removing bright galaxies 
and subtracting zodiacal and Galactic foregrounds from the measured sky brightness, we have 
successfully measured the CIB brightness and its fluctuations across a wide range of angular 
scales from arcminutes to degrees. The measured CIB brightness is consistent with previous 
results reported from COBE data but significantly higher than the lower limits at 70 and 160 
/tin obtained with the Spitzer satellite from the stacking analysis of 24-/im selected sources. The 
discrepancy with the Spitzer result is possibly due to a new galaxy population at high redshift 
obscured by hot dust. From power spectrum analysis at 90 /mi, three components are identified: 
shot noise due to individual galaxies; Galactic cirrus emission dominating at the largest angular 
scales of a few degrees; and an additional component at an intermediate angular scale of 10 — 30 
arcminutes, possibly due to galaxy clustering. The spectral shape of the clustering component at 
90 /mi is very similar to that at longer wavelengths as observed by Spitzer and BLAST. Moreover, 
the color of the fluctuations indicates that the clustering component is as red as Ultra-luminous 
infrared galaxies (ULIRGs) at high redshift, These galaxies are not likely to be the majority of 
the CIB emission at 90 /xm, but responsible for the clustering component. Our results provide 
new constraints on the evolution and clustering properties of distant infrared galaxies. 
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1. Introduction 

Since the Cosmic Infrared Background (CIB) 
was observed by the COBE satellite, it has been 
known that a large fraction of radiation en- 
ergy in the Unive r se is released in th e infrared 
(jPuget et al.lll996t iHauser et al.lll998l ). A com- 



pelling explanation for the CIB in the far-infrared 
is the thermal emission of interstellar dust asso- 
ciated with high-redshift galaxies, heated by the 
internal UV and optical radiation from young stars 
and active galactic nuclei (AGN), while the near- 
infrared background is d ominated by starlight 
from extragalactic sources (Hauscr & Dwekll200l[ 



Kashlinskv 11200a lLagache et alj|2005| ). Thus, the 
far-infrared CIB is a powerful probe of dust en- 
shrouded star-formation and AGN activity. Mea- 
suring the CIB in the far-infrared may constrain 
dust emission powere d by pre-galactic object s 
at very high redshift ( Coorav fc Yoshidal 20041) 
such as the first generat ion of stars (jBond et al. 
19861 : Santos et al.l 2002 ). which are expected to 
contribute to the near-infrared excess of CIB 



( Matsumoto et alJl2003 > 



In Figure Q] we summarize the results of the 
CIB measurements to date over the entire infrared 
and submillimeter wavelength range from previous 



space missions; COBE, IRTS, ISO (Hauscr et al 



1998; iFixsen et al 



Matsumoto et al 



199 
20051: 



Finkbeiner et al.. 

2009). 



Juvela et al 



2000; 



For 

comparison, we also show various foreground 
emission components; zodiacal light (sunlight 
scattered by interplanetary dust), zodiacal emis- 
sion (thermal emission from interplanetary dust), 
Galactic light (unresolved stars), diffuse Galactic 
light (starlight scattered by interstellar dust) and 
Galactic cirrus emission (thermal emission from 
interstellar dust), and the CMB. Note at mid- 
infrared wavelengths it is currently impossible to 
detect the CIB because the zodiacal emission is too 
bright. In the near-infrared and far-infrared, the 
foreground emission is relatively weak, and care- 
ful modeling and subtraction of the foreground 
enables one to extract the CIB from the measured 
sky brightness. 

As seen in Figure [TJ the CIB spectrum at wave- 
lengths longer than 200 /itm has been well con- 
strained with the FIRAS instrument on COBE 
(|Puget et al.lll996t IFixsen et al.lll998l ). However, 
results of the photometric measurements at wave- 



lengths shorter than 200 ^m with the DIRBE in- 
strument on COBE are divergengent in the mean 
levels of the CIB brightness, mainly due to the 
strong and uncertain foreground contamination of 
zodiacal emission, which dominates the sky bright- 
ness over all of the sky even though the Galactic 
foreground m ay be sufficiently weak in low cirrus 
regions (Hauser et al. Il998 : Lagache et al. 2000a : 



Finkbeiner et al.ll2000h . Although the CIB bright- 



ness was recently estimated using the ISOPHOT 
data, independently from the COBE data, the 90 
/zm data gave only an upper limit, and the mea- 
surement accuracy of the CIB brightness in the 
150-180 ixm range was in fact worse than that of 
COBE (j Juvela et al.ll2009h . Figure [Q clearly shows 
a wavelength desert of the CIB measurement at 
shorter far-infrared wavelengths, i.e., <200 pcax. 
Hence, new measurements of the mean level of CIB 
are demanded in this region. 

In the last decade, many observational efforts 
have been made to resolve the CIB into individ- 
ual galaxies by far-infrared deep surveys with 
infrared space telescopes su c h as I SO, Spitzer 
and AKARI dKawara et"aj] Il998l IPueet et af 



19991; iMatsuhara et allb oOO; iKawara et al. 2004 



Dole et all2004HFraver et al.ll2006al;lMatsuura et al 
120071 ; IShirahata et al.ll2008l) . and consequently the 
origin of CIB has become clear. As shown in Fig- 
ure Q] however, the detected galaxies account for 
only a small fraction (~ 10%) of t he measured 
CIB br ightness in the far-infrared. iFraver et al. 
(20063) claimed that they resolved more than 
half of the model CIB at 70/im into point sources 
in a single deep survey towards the GOODS-N 
field. In the mid-infrared, a lower limit of the CIB 
at 24 /urn is derived from the integrated number 
counts, and this is thought to account for ~ 70% of 



the model CIB a t 24 /im (jPapovich et al.l 12004) 



Dole et al.l (|2006l ) obtained lower limits for the 
CIB at 70 and 160 fj,m by stacking analysis of 24- 
fim detected sources, finding that the mid-infrared 
sources contribute ~ 80% of the CIB in the far- 
infrared, shown in Figure Q] by the dotted line. In 
the submillimeter range, similar stacking analy- 
sis of 24- /im galaxies against the deep surveys at 
250, 350 and 500 ^m by the Balloon-borne Large- 
Aperture Submillimeter Telescope (BLAST) ex- 
periment resulted in the conclusion that the 24- /im 
sources produce almost all the CIB in the submil- 
limeter measured with FIRAS ( Devlin et al. 2009; 
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Marsden et al. I l2009h . Although these studies us- 
ing the Spitzer 24-/mi surveys provided strong 
constraints on the mean CIB level, the current 
limit of direct measurement of CIB as diffuse emis- 
sion in the far-infrared range is still high enough 
to allow the existence of new population. 




wavelength [|im] 



Fig. 1. — Summary of the results of previ- 
ous CIB measurements in the infrared and sub- 
millimeter wavelength range with space mi s sions ; 



COBE/DIRBE ( diamo nds) (jHauser et al.l Il998t 
Finkbeiner etail l2000h . CO BE/FIRAS (shade d 



region in the submillimeter) (|Fixsen et al.l [1998). 
IRAS and ISO (thin line with downward ar- 
rows at 60 and 90 /jm and cross in 150-180 /xm) 
(|juvela et al .||2009|) , and IRTS (shad ed region m 



the near-infrared) (jMatsumoto et al.l 120051) . For 



comparison, we also show various foreground emis- 
sion components at dark sky; zodiacal light (ZL), 
zodiacal emission (ZE), Starlight (SL, K > 9 
mag), Diffuse Galactic light (DGL) and Galac- 
tic cirrus (ISD), and CMB. The integrated flux 
from the galaxy counts by deep surveys from the 
ground in the near-infrared and submillimeter and 
with the space telescopes in the mid- and fa r- 



1998: Pueet et al.l Il999t iMatsuhara et al.l 2000: 


Kawara et al.ll2004i iDole et al.l 


2004 iFraver et al. 


2006a; Matsuura et al.l 12007 


; Shirahata et al. 



2008), are indicated by the filled triangles con- 
nected with thin lines. Stacking of the 24 /im 
galaxies for the Spitzer/MIPS map (open trian- 
gles) and the BLAST map (open circles) results 
in good agreement with the predi cted CIB level 



from a galaxy evolution model by lLagache et al 
(|2004l ) (dotted line) 



Measuring the spatial fluctuations (anisotropy) 
of the CIB is a powerful method to investigate 
the unresolved galaxy population below the de- 
tection limit with little contamination from the 
foreground. The angular power spectrum of the 
CIB fluctuations is an important observable to 
trace the distribution of star- forming galaxies with 
respect to the clustering bias in dark matter halos. 
The fluctuation measurement is especially effec- 
tive at longer wavelengths, where direct measure- 
ment of the clustering of resolved galaxies is ham- 
pered by galaxy confusion due to the diffraction- 
limited resolution of current telescopes. In fact, 
power spectrum analysis of the CIB fluctuations 
has already been done for large area surveys in low 
cirrus regions; the FIRBACK survey at 170 
by ISO PHOT (|Puget et all Il999l: lLagache et af 



2000ri ). the MIPS / Spitzer data at 160 u rn in the 



agac he et al.ll20071) . and the 



SWIRE/GTO field _ . 

BLAST experiment (|Viero et al.ll2009l) . These re 



suits provided us with useful constraints on galaxy 
evolution and clustering. Information on the spec- 
tral energy distribution (SED) of the CIB fluctu- 
ations would be a further important clue in the 
investigation of the CIB properties. However, the 
quality of previous CIB fluctuation measurements 
at shorter far-infrared wavelengths, such as the 
EL AIS /Lockman survey at 90 ^m by ISOPHOT 
( Matsuhara et al.l 120001 ) and new analysis of the 



IRAS map at 60 ^m ([Miville-Deschenes et al 



2007), is not sufficient to add new information 
over the MIPS and ISOPHOT results at 160 and 
170 /im. 

In this work we aim to make the best mea- 
surements of the absolute brightness and fluctu- 
ations of the CIB since COBE, especially in the 
shorter far-infrared wavelength range which is as 
yet not extensively explored. This measurement 
is achieved using the high sensitivity of AKARI 
for diffuse far-infrared emission, high angular res- 
olution, and intensive selection of the survey field 
that provides the minimum foreground contami- 
nation. In the next two sections we describe the 
observational method and the data processing to 
produce the final calibrated image of the diffuse 
emission. In Section 4, we describe how the fore- 
ground subtraction so that the absolute brightness 
of CIB remains as an isotropic residual. In Section 
5, we describe the power spectrum analysis of the 
final image used to decompose the CIB fluctua- 
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tions into galactic and extragalactic components. 
In Section 6, we show the results of the data anal- 
ysis compared with previous studies of the CIB 
with scientific discussion. 

2. Observations 
2.1. Field selection 

The AKARI observations of the CIB were car- 
ried out for a far-infrared cosmological survey 
program (PI: S. Matsuura) as one of core scien- 
tific programs of the A KARI project, known a s 
Mission Programs (MP) (jMatsuhara et al.l 1200(1 . 
The survey field named the AKARI Deep Field 
South (ADF-S) is one of the lowest cirrus regions 



with contiguous area on the sky (jSchlegel et al 



1998) with an area of 12 square degrees near the 
South Ecliptic Pole; the field center is located at 
RA=4 /t 44 m 00 s , DEC=-53°20'00" (J2000). The 
survey field includes a ne arby cluster o f galaxies, 
DC 0428-53/ Abell S0463 (|Dresslerlll980h . The ob- 
servation log of this survey is summarized in Ta- 
ble ED 

The ADF-S field is one of the the best cosmolog- 
ical windows because of its low foreground emis- 
sion; not only a minimum of galactic foreground, 
but also low zodiacal background. High visibil- 
ity of this field for a polar-orbiting satellite such 
as AKARI is also important to achieve high sen- 
sitivity from long exposures. By comparison, the 
Lockman hole, the most popular low cirrus region 
as a cosmological window, has a relatively high zo- 
diacal brightness and low visibility due to its low 
ecliptic latitude. 

2.2. AOT and parameters 

In order to make a mosaic map of the entire 
field, a total of 200 slow-scan observations were 
obtained by combination of the Astronomical Ob- 
serv ation Templat e s (AOTs), FIS-01 an d FIS - 



02 (Ka wada et al.l l2007t iMatsuura et all |2007) 



These AOTs provide for multiple pointings as 
footprints of the field-of-view (FOV) along a scan 
path. The area covered by single observation of 
FIS-01 and FIS-02 is 12 arcminx0.6 deg and 8 
arcminxl.25 deg, respectively, for an exposure 
time of 10 minutes. Data in four photometric 
bands, centered at 65, 90, 140 and 160 /im, were 
simultaneously taken by a single observation, and 



the same sky is observed in all the bands except 
for a small difference in the FOV. 

The AOTs include the calibration sequence, 
i.e. dark current measurement by closing the cold 
shutter and responsivity check with an internal 
calibrator, during maneuvers before and after the 
observation and at the turning points of the round- 
trip scans. Such a highly redundant calibration 
sequence enabled us to correct any responsivity 
change, referring to astronomical calibration data 
taken by the same AOT. After the first maneuver 
to the target position, the cold shutter was opened, 
the sky signals are monitored during a settling 
time for the attitude control while the telescope's 
line of sight is oscillating around the target po- 
sition, and then the slow-scan observation starts. 
The signal monitoring before the observation is 
useful to check for any signal drift including time 
variation of any stray light (which is described in 
the next section). 

In order to produce a mosaiced image of the en- 
tire field with long-strip image segments, a total 
of 200 slow-scan observations with field overlaps 
of half a FOV in the cross-scan direction were car- 
ried out. The survey was completed in one and a 
half years over three observational seasons as sum- 
marized in Table [T] Mosaicking of the data taken 
in different seasons is not straightfoward, and re- 
quired careful treatment for the time variation of 
the zodiacal brightness. 

2.3. Supplemental observations 

In addition to the ADF-S survey, supplemen- 
tary observations using AOT FIS-01 were also 
made over small areas towards the Lockman Hole 
(LH) and the North Ecliptic Pole (NEP), as shown 
in Table [T] One of the LH observations was a 
mini survey to map over ~1 sq.deg for the perfor- 
mance v erification (PV) soon after the launch of 
AKARI (IMatsuura et al.l 120071 ). and another was 
included in the same MP program as the ADF-S 
to measure the galaxy confusion and to check the 
reproducibility of the measurement a year later. 
The NEP observation was a Director's Time pro- 
gram to evaluate the confusion noise in a relatively 
higher cirrus region compared to the ADF-S, the 
lowest cirrus region. As described in the latter 
sections, these observations are useful to check the 
validity of the foreground estimate and to evaluate 
the cirrus contribution to the sky brightness. 
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Fig. 2. — Sky brightness at 65 fim in ADF-S, 
Lockman hole and NEP fields, measured in vari- 
ous seasons (upper data points) and the residual 
brightness after subtraction of the zodiacal fore- 
ground (lower data points). The sky brightness 
values are only for those dates when the fields were 
observed, rather than the whole length of the mis- 
sion. The inset shows a magnified view of the data 
boxed with the dotted line. 
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3. Data processing and calibration 

In this section, we describe the methods of data 
reduction and analysis of the ADF-S data. The 
time series data are processed into a co-added 
image, using the FI S slow-scan analysis tool (SS- 
Tool Ver.20070914) (|AKARI FIS Data Users Manuall 
20071 ). The SS-Tool includes dark current subtrac- 
tion, absolute calibration, flat fielding, stray light 



removal, and co-addition of mult iple scans. The 
SS-To ol is described in detail in Matsuura et al.l 
(|2007h . therefore we describe them only briefly in 
the following sections. 

3.1. Time domain analysis 

The raw output of the FIS instrument is time 
series data, consisting of integration ramps for 
each detector pixel. Using the SS-Tool, we cor- 
rected the raw ramps for non-linearity and con- 
verted them to the average current values after 
removing spikes and glitches caused by cosmic-ray 
events and any other discontinuities. The SS-Tool 
process makes a dark current subtraction using the 
shutter-closed data before and after each observa- 
tion. 

The detector output changes with time depend- 
ing only on the sky brightness during the scan ob- 
servation. However, the output changes even when 
the telescope is pointed towards almost the same 
sky. This time varying component was identified 
as being from stray light caused by the Earth's 
limb emission illuminating the telescope baffle. 
The stray light appears in the case when the Earth 
avoidance angle of the telescope is smaller than 80 
degrees. The Earth avoidance angle changes with 
time depending on the orbital phase of the satel- 
lite and also with terrestrial season according to 
the inclination of the orbit even when the same 
field is being observed. The stray light appears 
typically for ~ 150 s at the begin and end of the 
observation time. At larger Earth avoidance an- 
gles stray light appears for ~ 200 s in the mid of 
the observation time. 

Stray light is removed in time domain by high- 
pass filtering the data with a window size of 90 
s, which was determined so to remove the slowest 
signal from the stray light but not to lose faster 
signals caused by astronomical emission. For the 
slow-scan observation, the time domain filter is 
equivalent to a spatial filter in the scan direction. 
We determined the window size of the filter in the 
scan direction as to be equivalent to the FOV size 
in the cross-scan direction, to ensure isotropy. 

3.2. Flat fielding 

Flat-fielding, i.e., correction of responsivity in- 
homogeneity in detector array, was carried out 
by measuring the mean brightness of the scanned 
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area for each observation. This method is based 
on the assumption that the sky brightness aver- 
aged over the observed area is the same for all 
array pixels. Accordingly, this method is effec- 
tive especially for dark sky dominated by zodiacal 
emission, whose spatial fluctuation is known to be 
very small from mid-infrared observations; their 
smoothness in arcmin scales being better than 1% 
((Abraham et al.lll998h . 

According to the results of previous far-infrared 
missions, e.g., IRAS and COBE, the sky bright- 
ness towards high galactic latitudes in the shorter 
AKARI bands at 65 and 90 /im should be dom- 
inated by zodiacal emission even at high ecliptic 
latitude; the expected fractions are 90% at 65/im 
and 70% at 90/zm, while at 140 and 160^m, Galac- 
tic cirrus emission is more likely the major (>50%) 
component in many cases. Since the smoothness 
of Galactic cirrus depends on the angular scale and 
mean brightness, the flat fielding in cirrus regions 
is not so excellent as the method using the zodia- 
cal emission. If we choose a sky with a brightness 
of about 10 MJy sr _1 , the smoothness of the cirrus 
can be better than 10%. 

In order to construct a flat-field template for the 
short wavelength bands, we have carried out ob- 
servations of zodiacal emission at ecliptic latitudes 
lower than 20 deg in the performance verification 
(PV) phase. For the long wavelength bands, we 
have also taken cirrus observations at intermedi- 
ate galactic latitudes of about 20 deg. However, 
the use of such a pre-measured flat has suffered 
from the responsivity change induced by radiation 
effects and the degradation of detector properties, 
which may not be perfectly corrected. 

For the ADF-S survey, in fact, flat fielding with 
the observed sky itself has yielded much better re- 
sults than the pre-measured flat, because the re- 
sponsivity change could be effectively corrected by 
the calibration sequence made during each obser- 
vation. Thus, we constructed such "self-flat" data 
by averaging the time series data for each pixel af- 
ter a 3-sigma rejection in time during the settling 
time for attitude control and the initial scan obser- 
vation. The flat-field accuracy with this method 
was generally better than 5%, and the actual effect 
of the flat field error for each observation was taken 
into account in the noise map derived from the sig- 
nal variance for multiple scans over the same sky. 



3.3. Absolute calibration 

3.3.1. Pre- flight measurement in laboratory 

Pre-flight measurements of the absolute gain in- 
cluding the detector responsivity were carried out 
in all the wave bands of AKARI in the laboratory, 
using an external blackbody source at various tem- 
peratures ranging from 17K to 30K with attenua- 
tors of 20-30dB. The measured gain factors show 
good agreement with each other to within ±5%, 
and we adopted a mean value as the final gain 
factor. 

In order to compare any measurements done at 
different times and different conditions either in 
orbit or in the laboratory, we have corrected the 
absolute gain value by referring to internal calibra- 
tor signals. The reproducibility of the calibrator 
signals was approximately ±5%, and this was also 
the limiting factor of the calibration accuracy. In 
the short wavelength bands at 60 and 90 fini, we 
estimated the accuracy of the pre-flight calibra- 
tion to better than ±7%, by combining the above 
uncertainties quadratically. 

While the monolithic Ge:Ga detectors for the 
short wavelength bands had fairly good uniformity 
of the responsivity over the array and showed a 
high stability in orbit, the stressed Ge:Ga detec- 
tor arrays for the long wavelength bands showed 
considerably larger responsivity differences among 
pixels and strong radiation effects in orbit. These 
properties made it difficult to correct for the re- 
sponsivity in the long wavelength bands even with 
the internal calibrator (c.f. flat field error) and 
caused additional calibration errors of ±6% and 
±9% at 140 and 160 /im, respectively. 

3.3.2. In-orbit measurement 

In order to measure the absolute gain of the 
detectors in orbit, we carried out observations of 
diffuse emission whose brightness had been pre- 
viously accurately measured with COBE, as de- 
scribed in the last section. 

For the short wavelength bands, we observed 
the zodiacal emission at low ecliptic latitudes. To 
derive the gain factor, we compared our observed 
data with the COBE/DIRBE model of zodiacal 
emission, which is reliable at low ecliptic lati 



tude to an accuracy better than 2% (jKelsall et al 
1998). In order to use this method avoiding 



6 



the contamination from Galactic and extragalac- 
tic emission, we differentiated the sky bright- 
ness measured at different ecliptic latitudes in low 
cirrus regions with similar HI column densities 
(Nffj ~2xl0 19 cm -2 ) and compared the differ- 
ential data with the zodiacal light model. The 
accuracy of this method is consequently limited 
by the calibration uncertainty of DIRBE ~10% 
(jHauser et alJll998h . 

For the long wavelength bands, we observed 
Galactic cirrus and compared the measured 
brightness with the COBE/DIRBE annual av- 
erage map (A AM). The sky brightness at low 
galactic latitudes in this wavelength range is dom- 
inated by Galactic cirrus, seasonal variation of 
zodiacal emission is smaller than 3%, and the 
AAM map is thus valid for this calibration. In 
addition to the calibration uncertainty of DIRBE, 
the low-frequency noise and instability of the de- 
tector during the observations increased the un- 
certainties to ±15% and ±19% at 140 and 160 
/.tm, respectively. 

3.3.3. Final gain factor and uncertainty 

For the short wavelength bands, the absolute 
gain factor of AKARI from pre-fhght measure- 
ments showed excellent agreement with that from 
the in-orbit calibration. Their difference was only 
5% at both 65 and 90 /im. It is noteworthy that 
this difference is smaller than either of the uncer- 
tainties of the pre-fhght calibration and the in- 
orbit calibration. For the long wavelength bands, 
the pre-fhght measurement gave a slightly differ- 
ent gain from that of the in-orbit calibration; the 
difference was 13% and 25% at 140 and 160 fim, re- 
spectively, but they agreed with each other within 
the uncertainties. 

In this work, we adopted our own gain factor 
from the pre-flight measurement to indicate the 
final result of the CIB brightness. For the fore- 
ground subtraction, however, we adopted the in- 
orbit calibration factor to compare our data di- 
rectly with the DIRBE result. In summary, the 
absolute calibration uncertainties of our CIB mea- 
surements taking all systematics into account are 
±7% at both 65 and 90 /zm, ±9% at 140 /mi, and 
±13% at 160 /jm. These numbers are sufficiently 
small enough for any significant detection of the 
CIB. More details about the absolute calibration 
are described in a separate forthcoming paper (S. 



Matsuura et al. 2010, in preparation). 

3.4. Co-added image 

Each slow-scan observation covers only a small 
patch of the sky corresponding to the scanned 
area. To produce the final mosaic image of the 
entire ADF-S field from the multiple scans, we co- 
added the time-series data in grid cells of the sky 
coordinates, according to the sky position of each 
detector pixel at each time of the data sampling. 
The co-addition process was done after subtract- 
ing the zodiacal foreground, which showed a sea- 
sonal variation as described in the next section. 

Even after subtraction of the zodiacal fore- 
ground, the mosaic image shows stripes as a scan- 
ning effect due to imperfect correction of the de- 
tector responsivity. In order to de-stripe the im- 
age, we determined the baseline levels of the zodi- 
subtracted images taken in the same season by a 
smooth filtering in the ecliptic longitude (cross- 
scan) direction with a filter width of 30 arcmin, 
which corresponds to three times the scan width. 
This filter works as a notch filter to reduce the 
noise only at the angular scale of the scan width 
preserving the fluctuation power at other angu- 
lar scales. The baseline levels are then subtracted 
from the zodi-subtracted time-series data. 

Finally, we produced the co-added image by cal- 
culating the weighted mean after a 3<r rejection 
in each grid cell of the image. We estimated the 
standard error of the mean value from the stan- 
dard deviation and number of data samples. The 
size of the grid cell is determined so as to con- 
tain sufficient number of data samples, according 
to the data interval, the scan speed, number of the 
scans, and number of the array pixels seeing the 
same sky. For the ADF-S survey, the grid cell size 
was set to be 20 arcsec x 20 arcsec each containing 
more than 8 data samples on average. 

4. Foreground subtraction 

In order to produce the final mosaic image of 
the entire field by the co-addition of the time se- 
ries data, we have adjusted the mean sky levels 
of multiple observations by subtracting the zodi- 
acal foreground with its seasonal variation. For 
accurate measurements of the CIB, estimates of 
the contribution from Galactic cirrus and exter- 
nal galaxies are also essential. In this section, we 
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describe our methods of the foreground subtrac- 
tion and show the final mosaic image of the entire 
ADF-S field. 

4.1. Zodiacal emission 

4- 1.1. Seasonal variation of zodiacal brightness 

Although the seasonal variation of the zodiacal 
brightness is sometimes an obstacle, it is rather 
useful to know the contribution of dust near the 
Earth to the zodiacal brightness. In fact, conven- 
tional models of the zodiacal cloud on the basis 
of the IRAS and COBE observations have been 
constructed relying on the local dust density and 
temperature derived from the s easonal variation 
(|ReacrJll98d iKelsall et al.lll998h . The amplitude 
of the seasonal variation is not affected by any 
uncertainty from the isotropic component of the 
zodiacal emission. 

In Figure [21 the measured sky brightness at 65 
/im at various locations in the ADF-S are plotted 
as a function of days after launch. Figure [3] is the 
same plot but at 90 /im. The sky brightness val- 
ues are only for those dates when the field were 
observed, rather than the whole length of the mis- 
sion. Statistical errors of the sky brightness from 
the RMS fluctuation of the image are of the or- 
der of a few kJy sr^ 1 (smaller than the size of the 
plot symbols). The sky brightness in the second 
season in Winter, about a year after launch, is 
significantly lower than that in the first and third 
seasons in Summer, as expected from a one-year 
period of the seasonal variation of zodiacal bright- 
ness. 

For comparison, the sky brightness towards the 
LH and NEP, which have been measured as sup- 
plemental observations in previously well explored 
fields, are also shown in Figure [5] and Figure [3] 
The LH for both PV and MP are brighter than 
the ADF-S because of the fields higher zodiacal 
brightness in spite of the similar density of Galac- 
tic cirrus. The NEP is brighter than the ADF-S 
due to the larger contribution from Galactic cirrus 
even though it has the lowest zodiacal brightness. 

4-1.2. Residual brightness after zodi- subtraction 

To measure the absolute brightness of the CIB, 
we subtracted the zodiacal foreground for each ob- 
servation by using the same COBE/DIRBE model 
as we used for the in-orbit calibration. The zodi- 
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55 50 45 40 E-lon [deg] 

Fig. 4.— Upper: Diffuse map of ADF-S at 90 
/zm after subtraction of the zodiacal foreground, 
drawn as linear contours. Brighter the contour, 
higher the surface brightness. The map size is ap- 
proximately 6 deg x 2 deg. Lower: The source 
removed map after masking out the sources down 
to ^20 mJy The masked area is filled with zeros 
to show the source locations clearly. 



acal brightness in the AKARI bands were calcu- 
lated by interpolation from the DIRBE bands with 
the color correction assuming a blackbody spec- 
trum. The zodi-subtracted data are composed of 
just the Galactic cirrus and CIB, as long as the 
zodi- model is correct. 

In Figure [H and Figure [3J residuals after sub- 
traction of the zodiacal foreground are compared 
with the raw data. The zodi-subtracted data for 
the ADF-S taken in different seasons show good 
agreement with each other, not only in the the 
mean levels but also in the detailed structure 
mainly due to Galactic cirrus. In addition, the 
residuals for the ADF-S and LH agree with each 
other except for the detailed structure, while the 
NEP is much brighter than the others, especially 
at 90 fim due to the larger contribution from the 
Galactic cirrus, as described later. At 65 /im the 
difference of the NEP brightness to the other fields 
is not so clear since the sky brightness is domi- 
nated by the zodiacal emission. The results indi- 
cate that the zodi-model is correct to an accuracy 
better than 5%, w hich is a typical v alue of the 
model uncertainty ( Kelsall et al. 1998f) . 
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4-1-3. Zodi- subtracted diffuse map 

The upper image of Figure |4] shows the final 
diffuse map at 90 fim produced by mosaicing the 
zodi-subtracted data. To obtain this mosaic im- 
age, destriping was made by first calculating the 
mean sky levels of the image strips taken by multi- 
ple scans after the removal of bright sources by a 3- 
sigma rejection method, followed by time-domain 
smooth filtering to the mean sky levels for the mul- 
tiple scans to correct any slight responsivity differ- 
ence that may not have been perfectly corrected by 
the calibration sequence. Finally all images are co- 
added. The window size used for the time-domain 
filter corresponds to three scan stripes in the space 
domain. A large contiguous area was successfully 
mapped, except for small lost data areas such as 
the inlets in the upper part of the map caused by 
problems with the ground station. As clearly seen 
in the map, there exists large scale structure from 
Galactic cirrus emission; with the central region 
darker than the surrounding area. Most of the 
bright spots can be attributed to nearby galaxies. 

4.2. Stars and galaxies 

Each star or galaxy in the field does not con- 
tribute much to the overall sky brightness, but the 
total contribution including faint sources down to 
the detection limit is non-negligible for accurate 
measurements of the CIB. Moreover, point sources 
seriously affect the sky fluctuations rather more 
than the absolute brightness. Hence, we masked 
out foreground stars and galaxies down to the 3a 
detection limit of ^100 mJy at 65 (im and ~20 
mJy at 90 fim, from the diffuse map. The mask is 
a circular area around each source with the radius 
normalized to be 1 arcmin at 100 mJy increasing 
proportionally with increasing flux. The lower im- 
age of Figure |4] is the resultant 90 /jm map after 
the source removal. The masked areas are filled 
with zero values to make the source locations clear. 

According to the point spread function (PSF) 
previously measured for stars and asteroi ds as 
calibration standards ( Shirahata et al. 2009), the 



remaining brightness in the out-of-mask region 
around a point source is less than 5% of the peak 
brightness, and the integrated flux in the out-of- 
mask region is 20% of the total flux. The inte- 
grated flux of the out-of-mask region of all sources 
down to the 3<r detection limit is less than half 



of the total source flux. Consequently, the mean 
brightness of zodi-subtracted map after removing 
the sources in such way is ~ 3% lower than the 
original value including the sources, and the sky 
fluctuation power is also an order of magnitude 
lower. We have also confirmed that the fluctua- 
tion power has the same value even if we change 
the mask aperture by a factor of ±50%. This im- 
plies the contribution of the out-of-mask fringes of 
PSF to the total sky fluctuations is negligible. 




0.4 0.6 
SFDIOOum [MJy/sr] 



Fig. 5. — Correlation between the AKARI 65 fim 
map and the SFD 100 fim map. The dashed line 
is the best fit to the data. 
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Fig. 6. — Correlation between the AKARI 90 fim 
map and the SFD 100 fim map, same as Figure [5] 



4.3. Galactic cirrus 

To estimate the Galactic foreground, we stud- 
ied the correlation between the ADF-S map and 
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the SFD dust map (jSchlegel et al.lll998l) . which is 
widely used as a reliable map of Galactic cirrus. 
This correlation may be scattered by point sources 
since the angular resolution of AKARI is much 
higher than that of the SFD map based on the 
IRAS and COBE data. In fact, many bright galax- 
ies in the field are resolved by AKARI, but in the 
SFD map they are left unresolved and recognized 
as part of the cirrus structure instead. Thus, we 
calculated the correlation for the zodi-subtracted 
map after removing bright galaxies and reducing 
the map resolution to be the same as SFD. 

Figure [5] and Figure |6] show the correla- 
tion between the SFD map at 100 /im and 
the ADF-S map at 65 and 90 /jm, respectively. 
The 90 /xm map correlates significantly with 
the SFD map with a correlation coefficient of 
7'=0.43, and the linear coefficient from the best 
fit line, shown in Figure |6] by the dashed line, is 
I(90/im)/I(100/xm) = 0.7030±0.0007. This color 
ratio corresponds to a dust temperature of T=18.2 
K for n=2, where n is the emissivity index and 
I{v) ~ B(T) * v ~ n , and is consistent with typical 
dust properties (jSchlegel et al.lll998l ). The corre- 
lation for the 65 /xm map was marginally seen as 
r=0.05 with Z(65/im)/J(lQ0//m) = 0.283±0.002. 

It is remarkable that the zero intercepts of the 
best-fit lines for the correlations for both the 65 
and 90 /im data show significantly finite values, 
which implies a positive detection of an isotropic 
emission component of the CIB. To estimate the 
cirrus brightness in the long wavelength bands, we 
extrapolated the 90 /im brightness with a black- 
body spectrum of T=18.2 K for n=2, because the 
map quality was not high enough to allow the same 
correlation analysis as made for the short wave- 
length bands. 

5. Power spectrum analysis 
5.1. Method 

In order to carry out a fluctuation study of CIB, 
we have made a power spectrum analysis using 
the two-dimensional (2-D) FFT for the final zodi- 
subtracted maps of the ADF-S. For the 2-D FFT 
analysis, we first produced a co-added image of 
the entire field with Cartesian coordinates by a 
simple tangential projection. The image size is so 
small that the distortion by the projection does 
not affect the result (while a spherical harmonics 



expansion with multipoles is required for power 
spectrum analysis at angular scales larger than 10 
deg). Next, we subtracted the mean offset level 
from the image, we set the blank data area out of 
the boundary of the observed area (dark area in 
Figure HI to zero value, and we also set the source- 
masked regions to zero value. Then, we applied a 
2-D FFT to the entire image. Finally, we corrected 
the raw power spectrum in the spectral domain 
by correction factors such as effective area, noise 
spectrum and PSF, as described below. 

5.2. Effective area 

As seen in Figure 0] the survey field has a non- 
rectangular shape, and the image includes zero- 
value data. Such an irregular shape may distort 
the power spectrum, and the non-negligible num- 
ber of zero- value data points decreases the ampli- 
tude of the power spectrum. In order to obtain 
the real power spectrum avoiding such artifacts, 
we have corrected the power spectrum by multi- 
plying it by a correction factor corresponding to 
the effective area of the field. 

To derive this correction factor, we used a ran- 
dom noise simulation method by generating noise 
images of the same shape and area as the survey 
field, so that the root-mean-square noise of all the 
images have a constant value. Then, we calculated 
the power spectra of the noise images by using the 
2-D FFT. Since the simulated noise image should 
have a flat power spectrum, i.e., white noise, we 
can find any effects due to image irregularity as 
a change in distortion and amplitude of the noise 
power spectra. The correction factor was actually 
calculated as the ratio between the original white 
noise spectrum and the mean power spectrum of 
the simulated noise images. 

5.3. Angular power spectrum 

Figure [7] shows the power spectrum of the 90 
/im map including bright galaxies as a function of 
angular wavenumber. The amplitude is corrected 
by the effective area as described above. To ob- 
tain the one-dimensional power spectrum, we have 
projected the 2-D power spectrum onto a radial 
axis of polar coordinates in angular wavenumber 
by calculating the ensemble mean for phase angles 
in a coaxial ring of each bin of the radial axis. The 
bin width is approximately 0.003 arcmin -1 corre- 
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Fig. 7. — Raw power spectrum of the zodi- 
subtractcd ADF-S map at 90 fim including bright 
sources (filled triangles), and the PSF corrected 
power spectrum after subtraction of the instru- 
mental noise (thin line). The diamonds represent 
the measured power spectrum of the instrumental 
noise, and the dotted line is the result of a poly- 
nomial fit to the noise spectrum. 



sponding to the inverse of the angular size of the 
map. The error bars of the power spectrum denote 
the standard error calculated from the standard 
deviation and number of data in each radial bin. 

5.4. Noise estimate 

In order to obtain the fluctuation power with 
astronomical origin, the subtraction of the instru- 
mental noise is necessary. We estimated the noise 
power spectrum from the difference of two subset 
maps taken in the same field but in different sea- 
sons. The obtained noise power spectrum shown 
in Figure [7] was dominated by the low- frequency 
noise as fitted by polynomials with terms of 1// 
and l// 2 noise. For the noise subtraction we used 
the result of polynomial fit rather than the data 
points to avoid adding any artifacts. 

The estimated noise is mainly due to instru- 
mental noise but may include systematic errors of 
the zodi-model uncertainty and the responsivity 
correction. The seasonal variation of small-scale 
structures of the zodiacal emission may also con- 
tribute to the noise, but its contribution is also 
subtracted from the measured power spectrum as 
a noise component. This noise level gives an up- 
per limit for the residual fluctuations of zodiacal 



emission smaller than 0.5% of the mean level at 
angular scales of ~1 arcmin, while the previously 
repo rted upper limit at th e same angular scale is 
1% ([Abraham et al.lll998l ). 
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Fig. 8. — The normalized power spectrum of PSF 
at 90 una directly obtained from the ADF-S map 
itself by two methods. The thin line is the sta- 
tistical PSF obtained by differentiating the maps 
before and after removing point sources. The 
data points represent the average PSF for bright 
sources in the field. 



5.5. Beam correction 

The steep decrease of the fluctuation power 
towards large wave numbers is obviously due to 
low-pass filtering by the PSF, which has a beam 
size of approximately 30 arcsec in FWHM. In or- 
der to reproduce the intrinsic power spectrum of 
the map, we divided the raw power spectrum by 
the normalized power spectrum of the PSF con- 
structed from sources in the ADF-S field itself 
(self- measured PSF). Although the pre-measured 
PSF of bright stars as calibration standards could 
also be used for this purpose, observations of cal- 
ibration stars with the same AOT parameter as 
the ADF-S survey were not available because the 
calibration sources are too bright. 

Figure [8] shows two kinds of power spectra for 
the self-measured PSF; one is a statistical PSF 
spectrum from the difference between two spec- 
tra of the source-removed map and the original 
map, and the other is the average PSF spectrum 
of bright sources only in the field. To obtain the 
average PSF spectrum, we cut out postage stamps 
with a size of 40x40 arcmin 2 around the bright 
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sources with signal-to-noise ratios higher than 30. 
The statistical PSF spectrum is normalized at the 
small end of wavenumber, and the average PSF 
spectrum is scaled to fit the overall range. The 2 
methods show good agreement with each other in 
the range of 0.1-0.8 arcmin, though a large discrep- 
ancy appears at large wave numbers. Since this 
paper focuses on the power spectrum at large an- 
gular scales we have adopted the statistical PSF, 
because it actually provides a better result in or- 
der to reproduce the power spectrum of the map 
at large wavenumbers and it is also known that 
transient effects are st ronger for brighter sources 



(IShirahata et al.l l2009) 
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Fig. 9. — Power spectrum of the ADF-S map 
after removing point sources (filled diamonds) 
compared with the power spectrum including the 
sources (open circles). The dotted line is the third 
power law for Galactic cirrus. The dashed line 
represents the shot noise level due to individual 
contribution of unresolved galaxies. 



tions, as discussed in the following section. 

We have carried out a power spectrum analysis 
for the 65 fim data in the same way as the 90 /im 
data. However, the obtained power spectrum was 
not significant compared to the noise power, suf- 
fering from low sensitivity. It was also difficult to 
obtain a significant power spectrum at 140 and 160 
/im because of the presence of low-frequency noise 
originating from the responsivity change induced 
by cosmic-ray hits. 

5.6.2. Bright source removal 

The raw power spectrum in Figure [JJ is obvi- 
ously contaminated by shot noise with a flat spec- 
trum at larger wavenumbers, originating from ran- 
domly distributed bright galaxies in the field. To 
investigate the fluctuations originating from un- 
derlying faint galaxies (the building blocks of the 
CIB), we derived the power spectrum removing 
the point sources down to the detection limit of 
AKARI, -20 mJy. 

Figure [§] shows the PSF-corrected power spec- 
tra before and after the point source removal. 
The shot noise power, in contrast to the cirrus 
component, decreases to 10% of the value before 
the source removal. The increased power at — 1 
arcmin -1 may be the signature of small-scale cor- 
relation between galaxies and source confusion or 
may be due to the uncertainty of the PSF as seen 
in Figure [8] Study of the power spectrum at these 
small angular scales is reserved for future work. 
At the other wavelengths, we have only obtained 
upper limits for the shot noise, which are not so 
meaningful as new limits. 



6. Results and Discussion 



5.6. Resultant power spectrum 

5.6.1. Raw power spectrum 

The power spectrum after the noise subtrac- 
tion and the PSF correction is shown in Figure [JJ 
There exist two major components in the power 
spectrum. One is the constant component domi- 
nating at wavenumbers larger than 0.1 arcmin , 
and another shows a steep rise towards the small- 
est wave number following a fc 3 -power law. The 
former is regarded as the shot noise (Poisson noise) 
mainly due to bright galaxies in the field, and the 
latter can be identified as Galactic cirrus fluctua- 



6.1. Mean CIB levels 

6.1.1. Absolute brightness towards the ADF-S 

The resultant CIB brightness was measured as 
the zero-intercept of the correlation between the 
zodi-subtracted map and the SFD map as shown 
in Figure [10] and summarized in Table [2] with sta- 
tistical and calibration errors with a n uncertainty 



in the zodiacal emission model of 5% (Kelsall et al 



1998). In Table [2] the mean brightness of zodiacal 
emission and Galactic cirrus emission over the en- 
tire ADF-S field are also shown as components of 
the measured total sky brightness. The measured 
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brightnesses in Table [2] are all color-corrected as- 
suming a flat spectrum. Although the zero-point 
offset of the SFD map may cause a systematic un- 
certainty in the CIB brightness, such uncertainty 
has been estimat ed to be smaller than 0.1 MJy 
sr" 1 at 100 /im dFinkbeiner et al.l l2000h . which 



does not significantly alter the results. In conclu- 
sion we can say we have significantly detected the 
CIB in the short wavelength bands. We have also 
tentatively detected the CIB at 140 and 160 by 
extrapolating the cirrus brightness from 90 ^m to 
each wavelength assuming a T=18.2 K blackbody 
spectrum. 

In order to obtain a more robust result, if we 
adopt another zodi-model with a very powerful as- 
sumption about the celestial signal called "very- 
strong no-zodi" principle, i.e., no extrasolar emis- 
sio n at 25 /xm exists at high ecliptic latitudes 



(Wri ght et al.l 119981 : lWrightll200tt iLevenson et al 
20071 ). then the zodi-model brightness becomes 



~10% higher than what was assumed for our fore- 
ground subtraction (L. Levenson 2009, private 
communication). On this assumption, the residual 
background decreases completely to zero at 65 /im 
but still remains at 90 ^m at a significant positive 
value of 0.50 MJy sr _1 , because of the bluer color 
of zodiacal emission. This lower limit of the CIB 
at 90 /jm is shown in Figure 1101 In the long wave- 
length bands, the zodi-model difference changes 
the residual background by at most 10%, which is 
smaller than the statistical error. 

6.1.2. Isotropy 

In Table [H we summarize the zodi-subtracted 
sky brightness in various fields as a mean value 
with the error estimated from the brightness vari- 
ance for a number of observations in each field. 
The error does not include either the calibration 
error or the zodi-model uncertainty. The table 
also indicates the residual brightness after sub- 
tracting the cirrus brightness, which were esti- 
mated from the HI 20 cm line map by the Leiden- 
Argentine-Bonn (LAB) survey ( Kalberla et al.l 
2005). The cirrus brightness was obtained by 
assuming a conversion factor from the Hl-column 
density to the 100 /im brightness of, 0.6 ±0.1 MJy 
sr _1 /10 20 cm~ 2 , as the median value of various 
COBE-based results with their variation as an 
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Fig. 10. — The CIB brightness measured in the 
ADF-S (filled circles), compared with the previ- 
ous CIB limits from COBE/DIRBE (diamonds), 
COBE/FIRAS (shaded region), IRAS (open cir- 
cle with downward arrow at 60 /mi), ISO (thin line 
with downward arrow at 90 /im and cross in 150- 
180 /im), Spitzer (upward arrows) and BLAST 
(open circles at 250, 350 and 500 /im). The filled 
triangle at 90 /mi is the lower limit of the CIB de- 
rive d for the case of "very-strong no-z odi" princi- 
ple (|Wrightll200l1 : ILevenson et al.ll2007h . The solid 
line shows a typical spectrum of the zodiacal fore- 
ground at the ecliptic poles. The dashed line is 



the C MB. Model predictions from lLagache et al 



uncertainty ( Boulang er et al.1 11996c iReach et al 



1998; ISchlegel et al.l 119981: lLagache et a l. 2000a) 



(|2004l) and iTakeuchi fc Ishiil (|2010l) are indicated 
by the dotted and dot-dashed curves, respectively. 
The open triangles and the open squares are lower 
limits obtained from the integrated flux of the 
galaxy counts with Spitzer and AKARI, respec- 
tively. 



For the wavelength conversion from 100 /im to the 
AKARI bands, we used the correlation coefficients 
in Figure [5] and Figure [6] Since the angular resolu- 
tion of the HI data is 10 arcmin, much larger than 
that of our data, we used the mean brightness of 
our data and of the HI data in each field. The 
residual brightness after subtraction of the cirrus 
brightness in different fields agree with each other 
within the error bars at both 65 and 90 /im. 

The ratio of the Hl-column density to the far- 
infrared brightness assumed above may be an over- 
estimate, because it tends to be determined by 
the data in the brighter cirrus regions. In low cir- 
rus regions the ratio could be smaller due to the 
contribution of uncorrelatcd components, such as 
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dust associated with ionized gas. In fact, the mean 
residuals for the three fields in Table are slightly 
smaller than the results in Table [2] from the corre- 
lation of the ADF-S and SFD maps, though they 
are consistent with each other within the total er- 
ror including the calibration error. The contribu- 
tion of cirrus could be more accurately estimated 
in the long wavelength bands, where the cirrus 
emission is brighter. However, we have not made 
the same analysis as the short wavelength bands, 
because accurate measurement of CIB was possi- 
ble only in the ADF-S field which has highly sta- 
tistical samples. 

6.1.3. Comparison with previous observations 

As seen in Figure [TUl our results in the four pho- 
tometric bands show good agreement with previ- 
ous CIB measurements from COBE at 60, 100 and 
140 /ira ( Hauser et al. 19981) . They are also con- 



sistent with an upper limit fr om IRAS at 60 /Ltm 
(|Miville-Deschenes et aDl2002h and with an upper 
limit at 90 fim and the measured CI B level in the 



150-1 80 /xm range measured by ISO (jJuvela et al 
20091) . 



The integrated galaxy counts (iDole et al.l 2004 ; 



Matsuura et al.l 120071 IShirahata et al.l 120081 ) ac- 



count for only a few per cent of the CIB brightness. 
Stacking analysi s of the Spitzer 24 /urn sources 
( Dole et aLll2006l) substantially extended the lower 
limits at 70 and 160 /im to the pred iction from the 



galaxy evolution model presented in lLagache et al 
(|2004l ) . Our measured CIB brightness at 90 /im is 



obviously higher than the Spitzer result, and our 
result implies the existence of large number of un- 
known sources, which may have been missed in 
the stacking process since they were too faint to 
be detected at 24 /jm with Spitzer. 

A hypothetical source with such an excess 
brightness should have peculiar spectra with a 
peak around 100 /Ltm, so that they have only a 
small contribution to the CIB in the other wave 
bands. High levels of the CIB at 100 ixm have al- 



ready been seen with COBE (jHauser et al. 19981 : 



Finkbeiner et al.l 120001 ) but have been thought to 
be due to calibration errors and foreground un- 
certainties ( Dole et al.l 120061 ). It is remarkable 



that this work corroborates with the COBE result 
using a different instrument. 



6.1.4- Possible origins of the excess emission 

If the excess emission at 90 /im originates from 
an isotropic component of the Galactic cirrus, 
which may not be accounted for in the SFD map, 
it should be tightly constrained by the previously 
measured Hl-column density unless the gas-to- 
dust ratio is extraordinary. If the excess bright- 
ness at 90 fim of ^0.4 MJy sr _1 is associated with 
neutral Hydrogen, the corresponding Hl-column 
density must be Nhi^ 1 x 10 20 cm~ 2 . This is 
much higher than the measured Hl-column den- 
sity towards the ADF-S (see Table [3|) and unlikely, 
judging from uncertainty of the HI observation. 
Moreover, if the temperature of such isotropic 
dust emission is ~17 K, similar to the bulk of the 
Galactic cirrus emission, it should appear at longer 
wavelengths as a strong excess, but this is not sig- 
nificant in our data. Therefore, it is difficult to 
explain the excess brightness by dust emission as- 
sociated with neutral Hydrogen. 

Unexpectedly large amounts of ionized gas at 
high galactic latitudes as a warm ionized medium 
(WIM) is a another possible candidate for the ex- 
cess emission. In fact, dust emission from a WIM 
has been detected with the WHAM H Q survey, 
though the observed areas have been limited to 
relatively dense region s at low galactic latitudes 
(jLagache et al.1 l2000al ). Large differences in the 
residual brightness between different fields shown 
in Table |3] can be caused by the contribution of 
dust associated with a WIM. In fact, the ADF-S 
is darker than the LH in the far-infrared, while 
the Hl-column density towards ADF-S is higher 
than that towards LH. However, any WIM could 
be dust poor compared to the Galactic plane, be- 
cause grain destruction is expected to occur in low 
density regions of the interstellar medium, while 
the ionized Hydrogen column density to account 
for all of the excess brightness must be higher than 
that of the neutral Hydrogen. As in the case of the 
neutral Hydrogen, the dust temperature must be 
so high as to have an emission peak around 100 /im 
while not producing a significant excess at longer 
wavelengths, but the existence of such hot dust at 
high galactic latitudes has not been reported. Fu- 
ture measurements of the H Q emission and other 
interstellar lines in the ADF-S would be important 
tests of the WIM hypothesis to explain the excess. 

Star-forming galaxies at high redshift are an- 
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other possible origin of the CIB excess in the far- 
infrared. In Table [2j the measured CIB bright- 
ness is compared with the range predicted from 
a selection of galaxy evolution models, which 
take the contribution of a population of lumi- 
nous infrared galaxies peaking at z ~ 1 and 
ultra-luminous infrared galaxies at z > 2 into 
consi d eration (jLaeache et alJl2004t IPearson et al. 
20071 : iFranceschini et all 120081 : iTakeuchi fc Ishii 



2010! ) . Most of these galaxy evolution models 



give a lower brightness than our measured CIB 
at 65 and 90 /im, while the models agree with 
our result at 140 and 160 /im, the stacking anal- 
ysis of Spitzer data at 70 and 160 /jm, and the 
COBE results at wavelengths other than 100 /im, 
within the data uncertainti es. A model predic- 
tion from ITakeuchi fc Ishiil (|2010h in Figure [10] 
shows a fairly good agreement with the measured 
CIB at wavelengths shorter than 140 /im, but it 
gives too high brightness at longer wavelengths. 



especially in the submill imeter range ([Puget et al 



19961: iFixsen et al . 1998). It is obvious that our re- 



sult especially at 90 /im provides a tight constraint 
on galaxy evolution models. 

In contrast to the CIB measurement, the 
source counts at 90 /im with AKARI were 
much lower than those predicted from contempo- 
rary galaxy evolution models, while the galaxy 
counts at 65 and 140 /im arc consistent with 
the mode ls and the Spitzer re s ults at 70 an d 



160 /im (jMatsuura et al. 
20081 ). Such low counts at 



2007: Shirahata et al 



been found with ISO (Rodighic ro et al.l 2003; 



90 uui hay e also 



Heraudeau et al.l 12004 ). and it has been reported 



that many galaxy evolutio n models fail to repro- 
duce the 90 jj,m c ounts (jCharv fc Elbazj 12001 
ISerieant et al. I l2004l) . At 90 /im the integrated 
galaxy counts down to the detection limit of 
AKARI, ~20 mJy, can onl y account for less tha n 



10% of the CIB brightness (jMatsuura et al.ll2007l ) 
To reproduce the CIB brightness, a strong increase 
in the galaxy counts below the detection limit is 
necessary. If we assume a single power law with 
power index of -1.5 (Euclidian) for the cumulative 
galaxy counts, dN/dS ~ S' 1 - 5 , the power law 
needs to extend to 0.01 mJy. Such high counts re- 
quire an extraordinary evolution scenario because 
conventional galaxy evolution models predict a 
steep drop of the differential 90 /im counts to- 
wards the faint end below 1 mJy. 



To explain the previous COBE result of high 
CIB levels at 60 and 100 /tm, the possible con- 
tribution of a new population at high redshift 
(z — 2 — 3) with high dust temperatures (Td ~60 
K) heated by AGN, or accreti on onto a black 
hole, has already been discussed dFinkbeiner et al ' 



20001 iBlain fc Phillips! l20oi iBlain et al.ll2004l ). If 



this is true, our finding of the CIB excess at 90 
/im could also be explained by such high-z hot 
sources. Many results from the Spitzer 24 /im and 
X-ray surveys suggest that LIRG and ULIRG at 
z ~ 1 rather than AGN at high redshift consti- 
tute main bulk of the CIB, and the contribution 
of AGN to the CIB i s thought to be at most 10% 
(|Treister et al.l 120061 iDevlin et all 120091) . There- 
fore, the CIB excess is more likely to be com- 
posed of high- z hot sources of new galaxy pop- 
ulation rather than known obscured AGN, which 
could be detected at 24 /im. If the hot source is so 
heavily obscured that the dust extinction is strong 
even in the mid-infrared, and if the dust grains are 
so small that the emissivity index is a maximum, 
A~ 2 , at longer wavelengths, they can contribute 
to the CIB selectively at short far-infrared wave- 
lengths. 

6.1.5. Energetics and cosmological implications 

The total integrated energy of the CIB excess 
in a range of the AKARI bands is estimated to 
~10 nWm 2 sr _1 , which is comparable to the CIB 
energy asso ciated with the 24/iin galaxies of 24 
nWm 2 sr _1 ( Dole et al. 20061 ) . If we assume a con- 
stant star formation rate (SFR) over cosmic his- 
tory of z > 1, the excess energy corresponds to 
the SFR of - 0.1 Meyr^Mpc -3 , which is an or- 
der of magnitude highe r than that in local universe 
(jHauser fc Dwekl 120011 ). If the radiation energy is 
produced by a star burst at z, the total hydrogen 
mass fraction converted to metals in this epoch is 
estimated to AX ~ 0.004(1 + z). For redshift of 
z > 1, more than 1% of baryon mass density needs 
to be converted to metals. Starburst at high-z as 
the origin of the CIB excess may cause over pro- 
duction of metals beyond the metal abundance in 
local universe. Black hole accretion may pray im- 
portant roll to g enerate the background radiation 
at high redshift (jFinkbeiner et alj|2000h . 

In the near-infrared, the CIB excess which can- 
not be explained by the integrated flux of known 
galaxy populations has been reported, and the ori- 
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gin of the excess as first stars, or proto-galaxies, 
at z ~ 10 has been discussed in terms of the spec- 
tral signature of the redshifted Ly-q in the CIB 
near 1 /mi, as seen in F igure HI (iKashlinskv Il2005t 
Matsumoto et al.l 120051 ) . Findings of the CIB ex- 
cess in both near- and far-infrared lead us to con- 
sider their connection; the UV light from the first 
massive stars which produce the near-infrared ex- 
cess might be partly absorbed by the self-produced 
dust, re-emitted in the mid-infrared and redshifted 
into the far- infrared band at present. Accretion to 
the first blackhole as a remnant of the first star 
explosion might also be the hot source at high 
redshift. Dust production by the first stars and 
successive process as blackhole formation and ac- 
cretion is little understood, and this research sub- 
ject would be proceeded by future observations of 
infrared galaxies and the CIB in many details. 

6.2. Fluctuations of the CIB 

6.2.1. Galactic cirrus 

In order to obtain information on the CIB fluc- 
tuations from the measured power spectrum, it is 
necessary to estimate the fluctuation power of the 
foreground such as the zodiacal emission and the 
Galactic cirrus. Zodiacal emission, the most dom- 
inant foreground component, has been known to 
have a very smooth distribution, and its fluctua- 
tions in the mid-infrared have been measured with 
ISO to be smaller than 1% of t he mean bright 



tions can be expressed as 



p 

1 r/i 



ness at angular scales of 1 arcmin ([Abraham et al 
1998). This value is already much smaller than 



the shot noise level by unresolved galaxies. Al- 
though an unexpected time varying residual of 
zodiacal emission after the zodi-model subtrac- 
tion may contribute to the power spectrum over 
a broad range of angular scales, it has already 
been subtracted as a part of the noise power es- 
timated from the difference of the two sub-maps, 
as described in the last section. Thus, the main 
foreground is Galactic cirrus, especially at large 
angular scales. 

A A: 3 -power law of the fluctuation spectrum for 
Galactic ci rrus has been reported by m any au- 
thors, e.g.. IMiville-Deschenes et al. I d2007h . but it 
is remarkable that in this work the same power 
law was confirmed at such very low flux levels at 
90 /xm. The power spectrum of the cirrus fluctua- 



,(*) = P {k/ko) 



(1) 



where k is the angular wavenumber in arcmin , 
and Pq is the power normalized at ko = 0.01 
arcmin - . 7 is the power law index, ranging from 
2.5 to 3.1 and is known to be close t o 3 for the opti- 
cally t hin case at low cirrus regions (ISchlegel et al 



1998 1 : IJeong et all 120051: IMiville-Deschenes et al 



20071 ). As shown in Table IH our result gives the 
cirrus power at 90 /im of P = (5.5 ± 0.5) x 10 4 
Jy 2 sr _1 , and the power at 100 fjan is derived to 
be P = (1.1 ± 0.1) x 10 5 Jy 2 sr" 1 from the cor- 
relation coefficient between the ADF-S and SFD 
maps. 

At high values of k the power law index may 
not always be the same as the value derived at 
low k, where the cirrus component dominates the 
total fluctuation power. In previous work with the 
ISO observations, it has been proved that a single 
power law for the cirrus power spectrum is valid up 
to k ^0.3 arcmin , though this result is certain 
only in relatively bright regions. Since there is no 
strong reason to change the power law index at 
higher k, we assume a single power law over the 
entire wavenumber range for the estimate of the 
cirrus power. 

It has been reported that the fluctuation 
power of the cirrus emission scales with the 
mean brightness as a power law with index 
2 — 3 depending on the brightn ess. According 
to IMiville-Deschenes et al.l (J2007J) , the cirrus fluc- 
tuation power at 100 /im in the low brightness 
regime (B < 10 MJy sr _1 ) is given by 



P = 2.7 x 10 b B 2V Jy 2 sr~ 



(2) 



where B is the mean brightness at 100 fim. If this 
is the case, our measured power of the cirrus fluc- 
tuations corresponds to the mean brightness at 90 
/j,m of 0.14 MJy sr _1 , which is slightly lower than 
the value in Table O The relation between the 
fluctuation power and the mean brightness may 
not be valid under such low cirrus conditions. The 
result implies that it is dangerous to estimate the 
mean CIB brightness from the fluctuation power. 

6.2.2. Shot noise of unresolved galaxies 

We have derived the shot noise of unresolved 
galaxies from the source -removed power spec- 
trum in Figure [9] as the mean power level in a 
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range of 0.2-0.8 arcmin -1 , where the PSF is rela- 
tively accurate, to be P s hot = (3.6 ±0.2) x 10 2 Jy 2 
sr -1 at 90 /im. In Table |4] we summarize our re- 
sult for the foreground estimates at small angular 
scales, where the shot noise dominates, into known 
emission components. The fluctuations of zodiacal 
emission was estimated to be negligible from the 
discussion in the last section. The cirrus compo- 
nent following a & 3 -power law is negligible at such 
small angular scales (k > 0.2 arcmin -1 ). There- 
fore, the remaining shot noise can be attributed to 
the fluctuations from unresolved galaxies. 

In Table|4]we compare our result with the previ- 
ous w ork at 90 /im with ISOPHOT ([Matsuhara et al 
2000h and with a range of the shot noise lev- 
els predicted from various galaxy evolution mod- 
els. Since the detection limit of AKARI for point 
sources is better than that of ISO, the AKARI re- 
sult shows obviously lower fluctuation power than 
the ISO result. We have estimated the shot noise 
level for each galaxy evolution model by simply 
integrating the squared flux, S 2 , with the differ- 
ential source counts, dN/dS, as 



Pshot = / S 2 (dN/dS)dS 



(3) 



for fluxes below 20 mJy. The measured shot noise 
power is 3-5 times lower than any of the model 
predictions clearly due to the high source counts 
of the models near the detection limit. The result 
implies that the main bulk of the CIB at 90 /im 
consists of sources much fainter than the detection 
limit. 

Since the shot noise is the squared flux inte- 
grated with the source counts, it is more sensitive 
to the number density of bright sources than the 
mean CIB level, that is a linear integration of the 
flux. As described in the last section, the mean 
CIB level provides constraints on the source counts 
at the very faint end below the point-source detec- 
tion limit. In contrast, the shot noise can provide 
constraints on the source counts near the detection 
limit. 

The measured shot noise level can be repro- 
duced by the measured source counts extended 
down to 1 mJy with a Euclidian power law (a 
constant for the differential counts). In such a 
simple scenario, sources fainter than 1 mJy have 
little contribution to the fluctuations but become 
the main contributor to the mean brightness. In 



order to resolve the bulk of the CIB into point 
sources, future infrared telescope with the detec- 
tion limits much better than 1 mJy at ~ 100 /im 
would be required. 




Fig. 11. — Sub- maps cut out from the selected 
positions in the ADF-S map at 90 /im; left: center 
region with lowest brightness, mid: a cluster of 
galaxies (Abell s0463), and right: relatively high 
cirrus region. The image sizes for all maps are the 
same, lxl deg 2 . 




Wave number: k [ arcmin" ] 



Fig. 12. — Power spectra of the three maps in 
Figure [11] Circles, triangles and diamonds rep- 
resent data points for the central, the cluster of 
galaxies, and the high cirrus regions, respectively. 
These power spectra are raw data without noise 
subtraction and the PSF correction. 



6.2.3. Cosmic variance 

Since the ADF-S covers a large contiguous area 
of 12 sq.deg, the field contains various regions of 
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Wave number: k [ arcmin 1 ] 



Fig. 13. — Power spectra for the sub-maps same 
as in Figure [T2l after the point source removal. 



different source and cirrus densities. The power 
spectrum shown in Figure |H] gives the average 
properties of the entire field, but it can be biased 
due to strong features of specific regions in the 
map. 

In order to check for any cosmic variance, we 
have carried out power spectrum analysis for sub 
maps of lxl deg 2 area cut out from various po- 
sitions as shown in Figure [IT] One is the central 
region with the minimum cirrus level and source 
densities, the second region contains a r ich clus- 
ter of galaxies DC 0428-53 (Abell S0463) (jDressler 
1980), and the third is a relatively high cirrus re- 
gion with a moderate source density. As shown 
in Figure [T2j the power spectrum obtained for 
the cluster region is about an order of magnitude 
higher than the others. The central region shows 
the minimum fluctuation power as expected, and 
the high cirrus region has slightly higher power 
than the central region, though the point source 
contribution seems to be more prominant rather 
than the cirrus fluctuations. This result implies 
that the power spectrum is severely affected by 
the point source contribution and biased to high 
source density regions. 

In order to estimate the effect of point sources 
on the power spectrum in detail, we have made 
the same comparison of the three sub maps after 
removing point sources. The result is shown in 
Figure [13J All of the power spectra at the differ- 
ent sky positions show excellent agreement with 
each other at an accuracy better than 5% in the 
overall range of angular wavenumber. This result 



implies that the power spectrum after removing 
point sources is not affected by cosmic variance 
such as the clustering of nearby galaxies and can 
be regarded as a universal function for the CIB 
fluctuations. 

6.2.4-. Clustering component 

In the measured power spectrum in Figure |H1 a 
fluctuation component of ~ 10 3 Jy 2 sr _1 is seen 
in the range of 0.03-0.1 arcmin -1 . Since the k- 
dependence of this component is moderate and 
apparently less steep than the A: 3 -power law, it is 
not likely to be attributed to either shot noise of 
unresolved galaxies or Galactic cirrus. No excess 
fluctuations of Galactic cirrus on such a small scale 
have been found in the high cirrus density region, 
where the cirrus dominates the sky fluctuations. 
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Fig. 14. — The final power spectrum at 90 
/im (filled circles) , same as in Figure [9l is com- 
pa red with the Spitze r/MIPS 160 /.im result 
in lLagache et al.l (|2007f ) (open diamonds and a 
shaded region at k > 0.2 arcmin -1 ) and their 
model of galaxy clustering (dash-dotted curve), 
Galactic cirrus (dotted curve) and shot noise 
(dashed line). The Spitzer result is multiplied by 
a factor of 0.01 to adjust the vertical scale. The 
solid thin line is the best fit to the AKARI data 
with a multiple power-law model (see text). 



A plausible interpretation of the fluctuation 
power at intermediate angular scal e is the cluster 
i ng of s tar-forming ga l axies 



(|2003l ). lLagache et al.l (|20071 ) have reported their 



e.g. 



Perrotta et al 



finding of galaxy clustering in the power spec- 
trum at 160 /im with Spitzer after removing point 
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sources brighter than 200 mJy. In Figure [Ml we 
plot the Spitzer result from iLagache et ah (2007) 
multiplied by a scaling factor of 0.01 with our 
power spectrum. Their power spectrum shows an 
excess component in the 0.01-0.1 arcmin -1 range 
that cannot be explained by the shot noise of in- 
dividual galaxies and the power-law spectrum of 
Galactic cirrus, agreeing with our result. They 
compared a galaxy clustering model (the dot- 
dashded line in FigureJTU) with the measured spec- 
trum and directly derived the bias factor for struc- 
ture formation in a dark matter universe to be 
b ~ 2.4. 

As a result of the comparison of our power spec- 
trum with the Spitzer result, the two spectra in 
k < 0.1 arcmin -1 agree with each other surpris- 
ingly well not only in the overall shape of the 
power law component at very low frequencies by 
the cirrus and of the galaxy clustering component 
but also in details such as wiggles at k = 0.01 — 0.1 
arcmin -1 and a dip around k ~ 0.02 arcmin , 
which may show an oscillatory structure of the 
baryonic matter distribution. Such concordance 
between the results of two independent observa- 
tions with different instruments and in different 
fields leads us to the conclusion that the cluster- 
ing components in both AKARI and Spitzer power 
spectra have the same origin. 

It has been reported that galaxy clustering 
similar to the AKARI and Spitzer results was 
clearly found by the BLAST experim ent in the 
submillimeter range ( Viero et al. I l2009h . Since the 
BLAST bands are on the Rayleigh- Jeans side of 
the dust emission in the rest-frame galaxy, BLAST 
can be more sensitive to high redshift galaxies 
than those of AKARI and Spitzer with less con- 
tamination by Galactic cirrus. The power spec- 
trum analysis for the BLAST map of 6 deg 2 area, 
found a strong signature of galaxy clustering with 
the power spectrum showing a fc 2 -power law over 
the range of k ~ 0.03 — 0.1 arcmin -1 in all the 
BLAST bands. Their power spectra show a very 
similar shape to those measured by AKARI and 
Spitzer, in terms of the power-law index of the 
clustering component and the angular scale at the 
knee between the shot noise and the clustering 
components, even though the limiting flux for the 
source removal for these three missions are dif- 
ferent; AKARI -20 mJy, Spitzer ~200 mJy, and 
BLAST -500 mJy. 



To extract the clustering component and to 
compare the AKARI result with the results of 
other missions on the same basis, we have also 
made a model fit to our power spectrum using a 
fc -2 -power law for the clustering component in ad- 
dition to the shot noise and the cirrus fluctuations, 
i.e., 



Pi{k)=P a 



P. 



hot 



P 



cluster i 



and 



Pcluster — Pcy^/^c) 



(4) 



(5) 



where fc c =0.03 arcmin -1 . Although such a sim- 
ple model cannot reproduce the complex structure 
seen in the measured spectrum as described above. 
For the best-fit case shown by the thin solid line in 
Figure [T4l the amplitude of clustering component 
is P c =(1.4±0.8)xl0 3 Jy 2 sr -1 . 

The fluctuation amplitude, 5I Vl of the cluster- 
ing component relative to the mean CIB level, 
< I v >, can be estimated from the power spec- 



trum at a given wavenumber to be (jKashlinskv 
20051) : 

SIu = V(k 2 P2(k)/2ir) (6) 

Using this relation, our measured fluctuation 
power of the galaxy clustering P c corresponds to 
a relative fluctuation amplitude of 8I v /< I v >= 
(2.3 ± 1.7) x 10 -3 . 

It is notable that the CIB distribution at 90 
/im after removing the point sources to a faint flux 
level is unexpectedly smooth. It is also noteworthy 
that the AKARI data are of high enough quality 
to detect such very low fluctuation levels. 

6.2.5. SED of the CIB fluctuations 

A study of the SED of the CIB fluctuations 
should give constraints on the redshift distribu- 
tion and SED of the underlying galaxies and pro- 
vide important information on their origin. Fig- 
ure [15] summarizes the CIB fluctuations mea sured 
with AK ARI, Spitzer dLaeache et all [2007) and 



BLAST (jViero et all l2009h . for shot noise levels 



at k — 0.5 arcmin - and for galaxy clustering at 
k — 0.03 arcmin -1 , with the fluctuation amplitude 
<5/„ . For comparison, the measured mean CIB l ev- 
els and the model CIB by ( Lagache et al. 20041 ) as 
in Figure [15] are also shown. 

If the clustering components measured by 
AKARI, Spitzer and BLAST originate from the 
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the CIB in the submillimeter ( Devlin et al. 2009; 



Pascale et al. 20091 : Viero et al 



2009). This result 



200 300 400 500 

wavelength [uiii] 



Fig. 15. — The spectral energy distribution (SED) 
of the CIB fluctuations meas ured with AKARI at 
90 /mi, Spitzer at 160 /mi ( Lagache et al.l 120071) 



and B LAST at 250, 350 and 500 /im (jViero et al 
2009) for the shot noise levels at k = 0.5 



arcmin -1 (filled triangles) and the galaxy cluster- 
ing at k = 0.03 arcmin -1 (open triangles). The 
fluctuation power is converted to linear bright- 
ness scale as the fluctuation amplitude, 51(h) — 
\f[k 2 P2{k) i '27r) . The fluctuations are compared 
with the absolute CIB levels, /„, same as in Fig- 
ureHOl AKARI (filled circles), COBE/DIRBE (di- 
amonds), Spitz er and BLAST (open circles). The 
model CIB by lLagache et al. I (120041) is also indi- 
cated (dotted line) as a reference SED of the CIB 
composed of the 24 sources, which are stacked 
with the Spitzer and BLAST data to obtain the 
mean CIB levels. The dashed lines with data 
points of the filled squ ares and the open squares 
are SED of Arp220 from lKlaas et al.l (|200lh . which 
is redshifted to be at z = 2 and vertically scaled 
to fit the CIB fluctuations data for the shot noise 
(upper) and the galaxy clustering (lower). 



same underlying galaxies, the combined spectrum 
should be comparable to the SED of the most dom- 
inant galaxy population. In Figure [15] we compar 
the combined spectrum of the CIB fluctuations 



with the SED of Arp220 taken from (jKlaas et al 
20011 ). which is redshifted to z = 2 and scaled to 
fit the data. Although there is a non-negligible 
discrepancy between the two spectra, the overall 
colors are very similar. Therefore, the clustering 
component is likely to be formed by such luminous 
infrared galaxies at redshifts of z > 1, which are 
all resolved at 24 /im and account for ~ 100% of 



is consistent in that the color of the clustering 
component is very similar to the color of the CIB 
in the Spitzer 160 /mi band and the BLAST bands 
and also consistent with the mean CIB levels in 
those bands which agree well with the model CIB 
for the 24/im-selected galaxies. 

It is noteworthy that the shot noise as shown in 
Figure [15] has a very similar spectrum to that of 
the clustering and to Arp220, in spite of the fact 
that the shot noise depends strongly on the lim- 
iting flux of the point source removal, which are 
quite different among the missions. This result 
implies that the brightest underlying galaxies in 
the field are tracing the clustering structure, be- 
cause the CIB fluctuation is more sensitive to the 
brighter end of the source counts than the mean 
CIB level, as discussed in the previous section. 

Another remarkable result on the SED of the 
CIB fluctuations is that the relative fluctuation 
amplitude at 90 /im, 5Igo/< /go >, is quite small 
compared with those at the other wavelengths, as 
summarized in Table [5] This is caused by both 
very red color of the CIB fluctuations and the high 
CIB level at 90 /mi exceeding the integrated flux 
of the 24/im-selected galaxies, which can account 
for the mean CIB levels at longer wavelengths. 
According to this result, the CIB should have a 
spectrum peaking around 90 /imif it is not to con- 
tribute much at 160 /im and not to be detected at 
24 /im. The underlying galaxies building the CIB 
should also be individually faint and numerous to 
keep the mean CIB levels high while keeping the 
CIB fluctuations low. 

Although hot high-z sources with far-infrared 
luminosities of the order of 10 13 L P) discussed in 



Blain et alJ ([2004J) can considerably contribute to 



CIB at 90 /im, such very luminous and less numer- 
ous galaxies would largely contribute to the CIB 
fluctuations rather than the mean CIB brightness. 
According to the measured CIB properties as dis- 
cussed above, the bulk CIB may consist of hot 
high-z sources of lower luminosity but higher sur- 
face density. A simple way of reducing the CIB 
fluctuations without the SED change is to assume 
higher redshifts and higher temperatures by keep- 
ing Td/ (1 + z) to be a constant, and this could be 
the case from our result. 

To obtain a definite conclusion to the origin of 
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the CIB excess in the AKARI band, further in- 
vestigation is required. To prove the hot high-z 
source hypothesis, large a perture space telescop e 



missions, such as SPICA ( Nakagawa et al. 2008), 



are required to resolve the far-infrared CIB into 
point sources to very low flux levels (below 0.1 mJy 
( Swinyard fc Nakagawal2009T) ). Spectral measure- 
ment of the continuous SED of the CIB over a wide 
infrared range would also be a powerful discrimi- 
nator of the galaxy population. Such observations 
can be simply done with a small aperture tele- 
scopes on small satellites or sub-orbital missions 
such as rocket experiments. Future space mis- 
sions to obtain a direct measurement of the mid- 
infrared CIB mitigating the contamination by zo- 
diacal emission require out-of-ecliptic missions or 
deep space missions beyond the asteroid belt (as 
the source of zodiacal dust), which will be power- 
ful enough to probe into the un derlying galaxies 



which compose of the CIB excess (jMatsuura et al 
2002b . 



7. Conclusions 

As the result of a cosmological survey in the 
ADF-S field, we have successfully measured the 
absolute sky brightness in the four photometric 
bands of AKARI at 65, 90, 140 and 160 fim. 
AKARI's high sensitivity for diffuse emission, the 
field selection to minimize the foreground, and 
careful foreground estimates have enabled us to 
detect the CIB as an isotropic emission compo- 
nent. The high angular resolution of AKARI was 
also essential for removing the foreground galaxies 
as point sources. 

The measured CIB brightness in all AKARI 
bands is consistent with the previous COBE/DIRBE 
results, and the 90 /im brightness is about twice 
higher than a simple interpolation from the Spitzer 
results at 70 and 160 //m obtained by stacking 
analysis of 24 /im-sclcctcd sources. The result 
suggests the exstence of new population(s) with 
high temperature spectra in addition to the pre- 
viously explored infrared galaxy population. 

By a power spectrum analysis of the 90 /jm map 
after removing the point sources down to ~20 mJy, 
we found the shot noise component due to the un- 
derlying galaxies. We also detected the signature 
of galaxy clustering at angular scales of a few tens 
arcmin. As a result of a comparison of our mea- 



sured power spectrum with the Spitzer result at 
160 fj,m and the BLAST result in the submillime- 
ter range on the same angular scales, the SED of 
the clustering component is found to be as red as 
that of ULIRGs at high redshifts, and the main 
bulk of the CIB may be composed of another pop- 
ulation with little contribution to the clustering 
power. 

Our results provide new information on the evo- 
lution and the clustering properties of infrared 
galaxies. To obtain a definite conclusion on the 
origin of CIB, further investigation is required, 
e.g., the study of the SEDs of infrared galaxies 
in detail and a cross-correlation study with multi- 
band images. Resolving the excess brightness into 
point sources and measuring the CIB spectrum 
continuously over a wide range would be the main 
target of on-going and future infrared missions; 
Herschel, SPICA and other satellites, deep space 
missions, and sub-orbital programs. 

We thank to L. Levenson for his help to cal- 
culate the zodiacal for eground in the far- infared 
based on the model in IWright et al.l ([19981 ). This 
research is based on observations with AKARI, a 
JAXA project with the participation of ESA. This 
work was supported by KAKENHI (19540250 and 
21111004). 
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Table 1: Observation log. 



Field 


Center coord. (J2000) 


Observation date 


AOT (parameter) 


N obs 


ADF-S 


4 h 44 m 00 s , -53°20'00" 


2006 Jul - 2006 Aug 


FIS-02(2.0;15;1) 


53 


ADF-S 


4 /l 44 m 00 s , -53°20'00" 


2007 Jan - 2007 Feb 


FIS-02(2.0;15;1) 


63 


ADF-S 


4 /l 44 m 00 s , -53°20'00" 


2007 Jul - 2007 Aug 


FIS-01(1.0;15;1) 


84 


LH-PV a 


10 h 50 m 36 s , 57°27'09" 


2006 May 


FIS-02(2.0:15;1) 


11 


LH-MP b 


10 h 50 m 25 s , 57°36'37" 


2007 May 


FIS-01(2.0;08;1) 


3 


NEP-DT C 


17 h 55 m 24 s , 66°37'32" 


2007 Oct 


FIS-01(2.0;08;1) 


4 



"Performance verification (PV) phase observations 
6 Mission program (MP) observations 
c Dircctor time (DT) observations 
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Table 2: Summary of foregrounds and mean CIB 
levels in ADF-S 



Wavelength 


Total a 


Zodiacal b 


Cirrus c 


CIB d 


Models 


w 


(MJy sr- 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


65 


4.31±0.04 


3.97 


0.09 


0.27±0.03±0.20 


0.1-0.2 


90 


3.57±0.02 


2.72 


0.22 


0.67±0.05±0.14 


0.2-0.5 


140 


2.36±0.16 


0.90 


0.57 


0.94±0.16±0.05 


0.5-1.4 


160 


2.02±0.20 


0.77 


0.60 


0.73±0.21±0.04 


0.5-1.6 



"Mean sky brightness of all observations derived with the in-flight calibration factor. The error denotes statistical error only. 

Mean zodiacal brightness of all observations calculated with the zodiacal emission model <Kelsall et al.lll998t) . 
c Mean cirrus brightness in the entire field converted from SFD map at 100/im ( Finkbcincr ct al. 2000) by using the linear 

correlation for the short wavelength bands and assuming the blackbody spectrum with T=18.2 K and n=2 for the long 

wavelength bands. 

''in the short wavelength bands, zero-intercept of linear correlation between AKARI and SFD maps to subtract cirrus brightness. 
In the long wave bands, the mean cirrus brightness was subtracted from the zodi-subtracted mean brightness. The gain factor 
measured in the laboratory is used. The first error includes both statistical error and calibration error, and the second error is 
uncertainty of the zodi-model. 

e CIB b r ightness range pr e dicted from galaxy evolution models, Lagachc ct al. (2004); Pearson ct al. (2007); Franccschin i et all 
ll200SD ; lTakeuchi fc IshiH (I2010T) . 
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Table 3: Field variance of CIB brightness 



Field 


N H i a 
(10 20 cm" 2 ) 


T b 

(MJy sr- 1 ) 


residual 


T b 
-'90 

(MJy sr" 1 ) 


residual 


ADF-S 


0.67 


0.329±0.003 


0.21±0.02 


0.835±0.004 


0.54±0.04 


LH-MP 


0.54 


0.311±0.057 


0.22±0.06 


0.873±0.021 


0.64±0.04 


NEP-DT 


3.73 


0.986±0.023 


0.33±0.09 


2.038±0.058 


0.42±0.22 


average 






0.22±0.02 




0.59±0.03 



a Taken from the LAB HI 20 cm survey l|Kalberla et al]|2005ft . 

''The mean brightness after subtraction of the zodiacal foreground with the errors estimated from the variance for a number of 
observations in each field. 

c Obtained by subtraction of the Hi-correlated component, assuming the conversion factor at 100 fim to be 0.6 ± 0.1 MJy 
sr - 1 /10 20 cm — 2 and extrapolating it to each wavelength by the blackbody spectrum with T=18.2 K and n=2. 
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Table 4: Small-scale fluctuations of foregrounds 
and the CIB a 



Zodiacal b 


Cirrus 


CIB P shot a 


Models ISO' 


< 0.6 x 10 2 


< 0.1 x 10 M 


3.6 ±0.2 x 10 2 


(1.1 - 7.2) x 10 3 2 x 10 4 



"Power pectral density at 90 /an in Jy 2 sr _1 . 

b < 1% of the mean brightness at angular scale of 1 arcmin (Abraham ct al. 1998). 
c Cirrus fluctuation power at k > 0.2 arcmin -1 estimated from the measured Pq. 
d The mean shot noise level in k = 0.2 — 0.8 arcmin -1 after point-source removal. 

"Galaxy evolution mo dels, lLagache et al. | i2004T) ; Ueong et all J2006fl ; iPearson et al] <2007fl : iFranceschini et al.1 <2008f) : 
iTakeuchi fc Ishiil j201Ch . 

/ISO result at 90 /im after removing galaxies brighter than 30 mjy. [Matsuhara et al. | J2000h . 
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Table 5: Measured SED of the clustering compo- 
nent of the CIB fluctuations 



Wavelength (fim) 


SIu (MJy sr- 1 ) 


5I u /< I v > 




Reference 


90 


(1.5±l.l)xl0 -3 


(2.3±1.7)xl0" 


-3 


This work 


160 


(1.8±0.2)xl0~ 2 


(3.2±0.3)xl0" 


-2 


1. 2 


250 


(1.7±0.3)xl0- 2 


(2.3±0.4)xl0" 


-2 


3,4 


350 


(1.3±0.3)xl0- 2 


(2.3±0.5)xl0" 


-2 


3, 4 


500 


(0.9±0.2)xl0" 2 


(2.3±0.6)xl0" 


-2 


3, 4 



Notes. The fluctuation amplitude of the clustering component estimated from Eq. 6 for k = 0.03 arcmin 1 . 
References. (l lDole et al.1 <2006ft : (2 lLagache et al.l j2007f ); O lMarsden et al.l teOOST ); (4 lViero et al.l j200Sft 
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